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Abstract 



We present the first results obtained with a Hybrid Molecular Dynamics algorithm 
applied to an N = 1 SU(2) Super- Yang-Mills on the lattice. We derive the Hamilton 
equations of motion for the system with Wilson gluinos and present preliminary results 
on small lattices. 



CERN-TH/96-122 
May 1996 



1 Introduction 



The supersymmetric extension of the Standard Model is now considered the most promising 
solution to the gauge hierarchy problem; this problem arises whenever the SM itself is consid- 
ered as the low-energy effective action of a more general theory including gravity. Moreover, 
it has recently been found that exact non-perturbative solutions of supersymmetric Yang- 
Mills theories can be found due to the peculiar properties of extended N = 2 supersymmetry. 
It is then interesting to begin the exploration of the non-perturbative aspects of supersym- 
metric theories by using their lattice formulation, in order to achieve a better understanding 
of the confinement and chiral symmetry breaking phenomena. The lattice formulation of 
a supersymmetric theory seems to be impossible since the discretization of the space-time 
explicitly breaks supersymmetry. Nevertheless, it has been shown ]l] that it is not necessary 
to put on the lattice an explicitly supersymmetric theory: much as in the case of chiral sym- 
metry, the main point is that the continuum limit of the theory should be supersymmetric. 
Thus, it is possible to simulate a non-supersymmetric theory with the same field content as 
the continuum one, with the requirement that in the a — > limit the supersymmetric theory 
is recovered. It has been shown (for N — 1 Yang-Mills in |Q and for N = 2 in ||) that 
the supersymmetric limit actually is strongly related to the chiral limit itself. In this way, 
in order to study N = 1 Super- Yang-Mills (SYM) it is only needed to introduce gluons and 
massive gluinos (fermion fields in the same representation as the gauge group of the gluons) 
and then recover the supersymmetric limit going to the m g = chiral one. It is evident 
from what was said above that in this theory it is fundamental to maintain an exact balance 
between bosonic and fermionic degrees of freedom. For this reason, an unquenched simu- 
lation should be performed, keeping also gluino loops. A recent paper [[| presented a first 
attempt to do an iV = 1 SU(2) Yang-Mills simulation on small lattices using the Luscher 
local bosonic algorithm. In order to perform a cross-check of the results, we have decided 
to adopt a different approach, namely the Hybrid Molecular Dynamics (HMD) algorithm of 

I- 

In section |2] we introduce the Wilson action for SYM theory and the Hamilton equations 
of motion; in section [3] we present our preliminary results, obtained on a small lattice, 
and finally in section |] we draw some conclusions discussing the feasibility of these kind of 
simulations and giving an outline of our future projects. 



2 Supersymmetric Wilson action on the lattice 

In this section the basic formalism, the notation and the main properties of the N = 1 
supersymmetric extension of the SU(2) Yang-Mills theory are introduced. The Hamilton 
equations for the Molecular Dynamics algorithm are also given. 

We recall that the parameter N refers to the number of ant i- commuting generators 
that are introduced in order to make a theory supersymmetric. The easiest supersymmetric 
extension of a non-Abelian gauge theory is iV = 1 SU(2) Yang-Mills, and so we will take this 
theory as our starting point for a study of the non-perturbative aspects of supersymmetric 
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theories with the lattice approach. The N = 1 Yang-Mills action contains two different 
pieces: the first is the usual Wilson action for the gauge bosons and the second represents 
the action term for the gauge fermions, the gluinos. We have adopted as fermion action 
the Wilson formulation in order to avoid the doubling of fermionic flavours. The full SYM 
action is 

SsYM — Sg + Sf, (1) 

where the standard Wilson action for the SU(2) gauge field is 

S 9 = ~ £Tr(C/ D ). (2) 
z □ 

The fermionic action should describe Majorana spinors in the adjoint representation of the 
gauge group. It is possible || to derive an expression formally identical to the usual Wilson 
action for the quark in QCD, replacing the link variables that are N c x N c matrices (where 
N c is the rank of the gauge group) with new fields in the adjoint representation V^, using 
(JVg — 1) x (jVj? — 1) matrices. We introduce the Majorana fermions as 

X(x) = T a X a (x) (3) 

with T a the gauge group generators, satisfying the trace condition 

Tr [T a T b ] = ^6 ab . (4) 

For SU(2) we can take T a = \o a , the Pauli matrices. We recall the gauge transformation 
properties of links and gluinos: 

Up(x) -> A(x)U f ,(x)tf(x + ft), 
X(x) -> A(x)A(x)A f (x). (5) 

The fermionic action, in the Wilson approach, is then given by 

Sf = ^£{A a (*)A a (x) 

[A q (x + v)V^(x)(l + lfl )X b (x) + X a (x)V^ b (x)(l - lfl )X b (x + //)]}, (6) 

where the adjoint link is given by 

V*\x) = l -Tz[Ul{x)a a U,{x)a b \. (7) 

It is easy to verify that belongs to 0(3), and so = Vj. In the continuum limit, the 
fermionic action @ gives the interaction Lagrangian 

Cj = i/ abc A a 4 7 ,A c . (8) 



The action can be rewritten in the following way: 



S f = -7 Tr 
3 4 



In(M^M)} , (9) 



where M is the fermionic matrix. The factor \ in eq. (|) contains a ~ due to the square root 
of det M^M and another | that takes into account the Majorana nature of the gluinos. A 
detailed description of how the gluino propagator and the n-point correlation functions can 
be obtained in the path integral formalism on the lattice can be found in ||. 

In order to simulate this theory, we have adopted the Hybrid Molecular Dynamics ap- 
proach introduced in [|J. Since the Super- Yang-Mills theory we want to reproduce contains 
just one gluino flavour, we are forced to use the so-called fi-algorithm. For this reason we 
cannot adopt the popular Hybrid Monte Carlo approach, which could be used only with a 
multiple of 4 Majorana flavours. As usual we introduce the Hermitian field P^ conjugate to 
the gauge field in order to write the Hamiltonian of the system as 

H = W^Pfa) + S g + S f . (10) 

The conjugate momenta P^ of the link variables are generated at the beginning of every 
Molecular Dynamics trajectory as complex Gaussian numbers distributed as exp(— |TrP^). 

Following |H we introduce a noisy estimator for (M^M) -1 in the form XX* i where x is 
the solution of the matrix equation 

M ] M X = <P (11) 

and = M^R. In this expression R is a vector of random Gaussian complex numbers with 
the same indices as the original gluino field, distributed as exp(— R*R). 

The Hamilton equations for the system are 

U^x) = tP^(x)U^x), 

P,{x) = l -(3T\G,{x)]- l -KT\F,{x)l (12) 

where the derivative is taken with respect to the Molecular Dynamics "time" and T is the 
projection operator on traceless anti-Hermitian matrices: 



T\B] = - 
[ 1 2 



(B - fit) _ _L T r(fi - fit) 



v c 



(13) 



The evolution of the conjugate momenta Pa fields is driven by the gauge "force 



Ul(x + fj, — v)U^{x — v)U v {x — u)] (14) 

and by the fermion "force" 

F,(x) = A^(x)Re[B b ;(x)}, (15) 
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where 



B b ;(x) 



a a Ul(x)a b U, 



x 



Tr 



(1 - ^)Y\x) X a *(x + //) + (! + 7 M )x fe (^)F-(x + /i) 



Dirac 



(16) 



and Y = M\- 

As in || , we assume that the sign of the Pfaffian operator (the square root of the fermion 
determinant) remains positive during the evolution of the system in the Molecular Dynamics 
"time" . In principle, we could take this phase into account by inserting it in the definition 
of the observables. 



3 Preliminary results 

In this section we present the first results we obtained with the Hybrid Molecular Dynamics 
algorithm for the Super- Yang-Mills SU(2) N — 1 theory. 

We have performed a test of the method using the same parameters as reported in || , 
namely /3 — 2.0 and K = 0.150. We have run the program on a lattice with the same size 
(4 3 x 8), obtaining compatible results for the plaquette, the absolute value of the Polyakov 
loop and the 2x2 Creutz ratio, defined in terms of k x I Wilson loops W(k, I) as 

_ H/(2,2)H/(1,1) 
X2 W(2,1)W(1,2)' { ' 

In all the runs the time-length of the Molecular Dynamics trajectories between different 
measurements was kept fixed at 0.5. This implies a number of leap-frog steps N step = 5, 10, 20 
for the three simulations performed. Our statistics is of the order of N m = 10000 measures 
for the full theory runs, and N m = 5000 in the pure gauge case. In this last case, we have 
used method 1 of p. With a short heat-bath run we obtain for the SU(2) pure gauge 
plaquette at (3 = 2.0 the value U PG = 0.5008(3). 

It is important to stress that the algorithm used introduces an intrinsic 0(At 2 ) error, 
where At is the Molecular Dynamics time step. It is then necessary to extrapolate to the 
At = limit in order to obtain physical results. Furthermore, the approach of the plaquette 
to the physical value is expected to be from above. Both behaviours, namely the At 2 
dependence and the approach from above can be seen in fig. [I] where the expectation value 
of the plaquette is presented function of At. 

We present in table [I] the results obtained for the three observables introduced above. 
The results are fully compatible with those reported in ||: for example the value of the 
plaquette extrapolated to At = is □ = 0.5058(4). We have no data for \P\ and \2 a t 
At = 0.025, but the values obtained at larger At are quite flat and already compatible with 

0- 

On the same observables, we have studied the autocorrelation times at different time 
steps, with the windowing method M. We build the connected autocorrelation function 
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At 


□ 




P 




X2 


N m 


0.100 
0.050 
0.025 


0.5098(3) 
0.5065(3) 
0.5063(3) 


0.0505(4) 
0.0500(4) 


0.5555(12) 
0.5559(13) 


10000 

8000 

10000 



Table 1: Average observables for f3 = 2.0, K = 0.150, V = 4 3 x 8 at different values of At. 
We show the plaquette, the absolute value of the Polyakov loop and the 2x2 Creutz ratio. 
N m is the number of measurements for each run. 



C(t) and from this an estimate for the integrated autocorrelation time 

When the ratio T/r int (T) is around 4 we should observe a plateau in the value of r int (T): this 
is exactly what happens, as can be seen in fig. ^| in a particular case. Results for integrated 
autocorrelation times are presented in table £2], together with a rough estimate of the errors. 



At 


T D 


r\ P \ 


T X2 


0.100 


7.5(5) 


1.0(2) 


2.8(3) 


0.050 


6.8(5) 


1.0(2) 


2.7(3) 


0.025 


6.7(5) 







Table 2: Autocorrelation times for the different observables at different values of At. 



Finally, we present a rough quantitative comparison of the relative speed of the two 
algorithms used for these simulations (namely, HMD and the Liischer local bosonic one 
adopted in 0); an extensive study of the scaling of the autocorrelation times with the 
different parameters relevant to the simulations has yet to be performed and, moreover, it 
is likely that further improvements can be implemented. Since the most time-consuming 
part of both algorithms is the application of the fermion matrix to a vector, we consider 
the "effort" E needed to obtain an independent configuration from a previous one as the 
number of times this application is performed, times the autocorrelation time of a particular 
observable (for definiteness, the plaquette). For each simulation, introducing N app as the 
number of times the fermion matrix is applied per step, we obtain 

E H aid = N step x N app xt. (19) 

For this preliminary study we have not yet tried to optimize the algorithm: in particular, so 
as to solve the matrix equation (|ll|), we used the Stabilized Bi-Conjugate Gradient without 
preconditioning. With this particular inverter and with the leap-frog scheme used to integrate 
the Hamilton equation, 

N app = 2 + AN iter , (20) 
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0.515 



Figure 1: Average plaquette versus At in the pure gauge case (squares) and within the full 
theory (circles). The upper cross at At = is reported by Montvay, while the lower one is 
the heat bath pure gauge result. The dotted line is the extrapolation to At = 0. 



,0000000 o 



I i , i L_ 

3 4 

T/T,n. 



Figure 2: Integrated autocorrelation time versus T/ri nt for the plaquette at At = 0.025. 



where Ni ter is the average number of iterations required to invert the fermion matrix. The 
convergence criterion we used for all our runs was \M^Mx — 4>\ 2 < 10~ 8 |x| 2 , obtaining N iter ~ 
25. For the most time-consuming simulation (At = 0.025, N ste p = 20), Ejjmd — 14000; for 
the whole set of three simulations at different At, we obtain Eh aid — 25000. 

For the local bosonic algorithm used in H], the "effort" is computed as 

E L B = N fidd X N app XT, (21) 

where Nfi e id is the number of bosonic field used in the simulation. Using the results reported 
in ||, we estimate (for the two most time-consuming simulations) Elb — 90000 for the one- 
step algorithm with 24 bosonic fields, and E^b — 25000 for the two-steps algorithm with 8 
fields [|. 

For a realistic study of the theory with the HMD algorithm, it is obvious that the inver- 
sion algorithm will have to be accelerated, for example by means of the standard red-black 
preconditioning, which we are currently implementing. Another important improvement 
of the algorithm can be achieved by using a higher-order Hamilton equations integration 
scheme, of the kind proposed in 0, in order to accelerate the convergence towards the 
At = limit. 



4 Conclusions 

In this paper, following the outline presented in recent papers |3], Q], we have begun an explo- 
ration of the feasibility of supersymmetric gauge theories simulations on the lattice. From a 
theoretical point of view, it is possible to study these kinds of theories on the lattice using the 
same approach adopted for chiral theories as usual QCD, namely the implementation of a 
non-chiral non-supersymmetric theory requiring the recovery of these symmetries in the con- 
tinuum limit. We have obtained results compatible with those presented by Montvay in || 
using a completely different algorithm, the Hybrid Molecular Dynamics one. From these first 
simulations it seems possible, with the present computer power and the present knowledge of 
fermionic algorithms (such as the local bosonic, the hybrid family or even with the "negative 
flavour number extrapolation" approach 0), to afford the study of non-perturbative aspects 
of supersymmetric extensions of gauge theories. 

The next step is to begin the investigation of the supersymmetric limit of the theory, 
which is expected to be reached jointly with the chiral limit PJ. In order to do that, we 
are planning to study the asymptotic behaviour of the scalar and the pseudoscalar gluino 
"meson". In the OZI approximation it is expected that the pseudoscalar mass approaches 
zero when going to the chiral limit, while the scalar should go to a non-zero value. 

Finally, we plan to investigate the behaviour of the gluino condensate. This quantity 
is expected to be a relevant parameter in the study of non-perturbative supersymmetry 

1 These numbers can be decreased to 44000 and 15000 respectively, by changing the mixture of heat-bath 
and over-relaxation sweeps in the update of the bosonic fields raj . 
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breaking: it is worth while to mention that a clear comprehension of the supersymmetry- 
breaking phenomena is of the uttermost importance if the supersymmetric extension of the 
Standard Model has to be a viable solution to the gauge hierarchy problem. In order to study 
this observable, it is necessary to understand its renormalization properties. It is perhaps 
important to introduce a non-perturbative renormalization scheme both on the action and 
on the operators [K|, which seems to be relevant in the study of the chiral (supersymmetric) 
limit of the theory with Wilson fermions. 
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